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ABSTRACT:  The  use  of  a  high-resolution,  ground-based  3D  laser  scanner  was  recently  evaluated  for 
terrestrial  site  characterization  of  variable-surface  minefield  sites  and  generation  of  surface  and  terrain 
models.  The  instrument  used  to  conduct  this  research  was  a  Leica  HDS3000  3D  laser  scanner.  The 
high-speed,  highly  accurate  ranging  system  has  a  360  deg  horizontal  *  270  deg  vertical  field  of  view  that 
delivers  positional,  range,  and  angular  (vertical  and  horizontal)  single  point  accuracies  (range  1  to  50  m) 
of  6  mm,  4  mm,  and  60  micro-radians,  respectively.  The  laser  is  a  class  3R  and  is  completely  eye-safe 
with  a  wavelength  of  523  nm  and  spot  size  of  <6  mm  at  a  distance  of  50  m.  The  pulse  rate  is  1,000  points/ 
sec  with  an  optimal  effective  range  up  to  100  m  which  is  capable  of  producing  a  maximum  point  cloud 
spacing  of  1.2  mm  in  the  horizontal  and  vertical  direction.  Two  study  sites  located  in  the  midwestem 
United  States  were  used  for  this  analysis.  A  very  dense  vegetation  site  (Grass  Site)  and  a  bare  soil  site 
with  intermittent  rocks  and  sparse  vegetation  (Dirt  Site)  were  selected  for  data  collection  to  simulate  both 
obscured  and  semi-obscured  minefield  sites.  High-density  scans  (range  0.2  to  2.0  cm  spacing)  were 
utilized  for  Cyra  target  acquisition  and  were  commensurate  with  size  and  distance  to  target  from  scanner 
location.  Medium-density  scans  (range  2.0  to  5.0  cm  spacing)  were  chosen  for  point  cloud  generation  of 
the  entire  site(s)  with  approximately  10  percent  edge  overlap  between  field  scans.  In  order  to  provide 
equivalent,  unobstructed  viewing  perspectives  from  all  scan  locations  at  each  site,  the  scanner  was  posi¬ 
tioned  on  a  trailer-mounted,  chain-driven  lift  and  raised  to  an  approximate  scan  height  of  7.6  m  above  the 
ground.  Final  registration  to  UTM  projected  coordinate  system  of  the  multiple  scan  locations  for  the  Dirt 
Site  and  Grass  Site  produced  mean  absolute  errors  of  0.014  and  0.017  m,  respectively.  The  laser  scanner 
adequately  characterized  surface  roughness  and  vegetation  height  to  produce  contour  and  terrain  models 
for  the  respective  site  locations.  The  detailed  scans  of  the  sites,  along  with  the  inherent  natural  vegetation 
characteristics  present  at  each  site,  provide  real-time  discrimination  of  minefield  site  components  under 
contrasting  land  surface  conditions. 
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Introduction 


Purpose 

This  measurement  and  analysis  effort  was  conducted  to  support  specific  key 
objectives  of  the  Engineer  Research  and  Development  Center’s  (ERDC) 
Countermine  Phenomenology  Program.  The  objectives  addressed  through  this 
research  primarily  include  the  development  of  an  improved  knowledge  of  geo- 
environmental  phenomenology  factors  impacting  both  mine  and  minefield 
detection  in  various  operational  environments  and  help  support  the  development 
of  algorithms  using  these  factors  to  improve  detection  capability  of  both  surface 
and  buried  mines  and  minefields. 

The  use  of  a  high-resolution  3D  ground-based  laser  scanner  was  evaluated 
and  assessed  for  characterizing  and  capturing  terrestrial  site  characteristics  of 
variable-surface  minefields  to  aid  in  surface  model  generation  and  ground  surface 
contour  mapping.  This  research  tool  is  one  of  many  being  implemented  by  the 
ERDC  Countermine  Phenomenology  Program  to  evaluate  potential  methods  for 
characterizing  and  delineating  background  features  within  minefields.  This  effort 
is  intended  to  support  and  improve  the  understanding  of  background 
phenomenology  with  respect  to  minefields  (Jackson  et  al.  2005). 


Background 

The  highly  accurate  and  dense  point  data  (or  point  clouds)  captured  by 
terrestrial  3D  laser  scanners,  such  as  the  Leica  HDS3000  system,  allow  for  the 
development  of  robust  datasets  for  GIS  modeling  efforts  and  dynamic  landscape 
visualization.  The  ground-based  instrument  is  closely  akin  to  its  airborne  light 
detection  and  ranging  (LIDAR)  counterpart.  However,  terrestrial  LIDAR  is 
acquired  with  more  efficiency  at  a  much  higher  resolution  from  a  more  stable 
platform.  Eliminating  the  need  to  correct  for  orientation  errors  common  with 
airborne  sensors,  terrestrial  3D  laser  scanners  produce  accuracy  measured  in 
millimeters  as  opposed  to  centimeters,  but  for  obvious  reasons  are  not  as 
effective  as  airborne  LIDAR  at  extracting  data  in  the  vertical  dimension,  such  as 
vegetation  height  and  ground  surface  elevation.  Other  limiting  factors  of  ground- 
based  LIDAR  are  its  restriction  to  small  geographic  areas  and  the  requirement  of 
numerous  scans,  which  can  become  labor-intensive  for  large  area  acquisitions. 

Multiple  laser  scanner  setups,  analogous  to  airborne  LIDAR  swaths,  are 
required  for  data  acquisition  when  implementing  terrestrial  LIDAR  and  must  be 
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coregistered  to  form  a  cohesive  point  cloud  representative  of  the  sampled  area. 
Laser  scanning,  in  general,  is  a  rapid  non-invasive  form  of  data  acquisition  that  is 
suitable  for  characterizing  areas  with  restricted  or  limited  access  or  where 
environmental  conditions  exist  that  may  limit  one’s  ability  to  gain  physical 
access  to  the  area. 

Airborne  laser  scanner  systems  are  abundant  on  the  market  today  and  the 
technology  is  considered  to  be  in  a  mature  state.  These  laser  systems  are  very 
complex,  being  more  a  geodetic  system  for  data  acquisition  and  more  a 
photogrammetric  system  for  data  processing  (Axelsson  1999).  What  is  lacking, 
however,  is  the  development  and  refinement  of  algorithms  and  data  interpretation 
methods  to  provide  various  surface  representations  of  the  scanned  area.  Likewise, 
little  is  known  about  the  implications  for  surface  modeling  from  the  use  of 
terrestrial  laser  systems. 

Terrestrial  LIDAR  —  being  active  sensing  devices  that  emit  electromagnetic 
energy  either  in  the  visible  or  near  infrared  part  of  the  spectrum  —  records  the 
amount  of  reflectance  and  laser  beam  return  time  from  the  scanned  feature.  The 
surveyor  can  define  the  area  of  interest  and  specify  the  angular  resolution  value, 
with  no  further  attendance  required  during  the  scanning  phase,  which  takes  only 
minutes  to  complete.  Linear  resolution  over  the  object  depends  on  distance, 
azimuth,  and  elevation  of  the  laser  beam,  as  well  as  terrain  surface  morphology. 
The  accuracy  of  point  positioning  is  a  function  of  distance,  the  number  of 
scanning  repetitions,  and  laser  spot-size  (Colombo  2003). 

Johansson  (2002)  documented  undesirable  artifacts  in  resulting  point  cloud 
data  generated  from  various  ground-based  laser  scanners.  Strange  effects  along 
edges  of  objects  and  difficulty  in  capturing  points  on  certain  surfaces  were  noted. 
The  author  emphasizes  small  spot  size,  good  range  capabilities,  and  high  posi¬ 
tional  accuracy  of  the  chosen  scanner  to  minimize  or  resolve  these  issues. 

Guamieri  et  al.  (2004)  described  the  insufficient  use  of  natural  objects  as 
control  points  for  model  georeferencing  and  emphasized  common  error  sources 
when  using  ground-based  LIDAR.  Detected  errors  were  directly  related  to  the 
scanner’s  accuracy,  intensity  response  from  the  reflected  laser  beam,  and  the 
operator’s  ability  to  identify  model  control  points.  The  authors  advise  the  use  of 
artificial  reference  targets1  whenever  possible  to  increase  the  accuracy  of  control 
point  selection.  These  fixed  targets  can  be  surveyed  and  will  therefore  lend 
proper  fitting  procedures  during  the  point  cloud  registration. 


1  Targets  defined  herein  as  scanner  reference  targets  (Cyra  targets)  used  for  control  point 
acquisition  during  laser  data  registration.  These  are  not  the  same  as  landmine  targets. 
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Processing 


Site  Description 

Two  variable-surface  minefield  sites  located  in  the  mid-western  United 
States  were  selected  as  study  site  locations  to  conduct  this  research  effort.  A  very 
densely  vegetated  site  (Grass  Site)  and  a  mostly  bare  soil  site1  (Dirt  Site),  with 
intermittent  rocks  and  sparse  vegetation,  allowed  for  characterization  efforts  to  be 
employed  for  obscured  and  semi-obscured  minefields.  Each  site  is  approximately 
rectangular  having  four-sided,  right-angled  polygons  with  opposite  sides  equal 
and  parallel  to  each  other.  The  dimensions  of  the  Grass  Site  and  Dirt  Site  are 
40  x  160  m  and  40  *  320  m,  respectively.  Both  minefield  sites  serve  as  test  beds 
for  deployed  mines  and  contain  M20,  Ml 9,  and  RAAM  buried  and  surface 
emplaced  mines  with  top  hat  and  Electro-Optical  Infrared  (EOIR)  red  board 
fiducial  markers  spaced  at  various  intervals. 

The  Grass  Site  has  about  a  5  to  10  percent  grade  with  a  northwest-facing 
slope.  It  is  comprised  mostly  of  thick  grass  with  varying  density  and  distribution 
over  the  field.  In  contrast,  the  Dirt  Site  is  relatively  flat,  with  the  exception  of  a 
drainage  ditch  oriented  north  and  south  across  the  center  of  the  field  on  the  east 
side.  The  Dirt  Site  has  been  plowed  several  times  and  consists  mainly  of  large 
boulders  and  smaller  rocks  with  sparse  patches  of  grass  vegetation.  The  basic 
capability  to  derive  contour  and  terrain  models  from  a  high-density  laser  data 
point  cloud  was  evaluated  for  these  two  geophysically  different  sites. 


Laser  Data  Collection 

Data  acquisition  took  place  during  late  July  and  early  August  2004.  A  Leica 
HDS3000  3D  laser  scanner  manufactured  by  Leica  Geosystems  HDS,  Inc. 
(formerly  Cyra  Technologies)  was  used  to  provide  a  high-definition  survey  of 
both  sites.  The  SmartScan  Technology™  of  this  unit  provides  a  maximum 
360  deg  horizontal  field  of  view  and  a  maximum  270  deg  vertical  field  of  view. 
The  scanner  emits  rapid  pulses  of  green  (523  nm)  laser  light  that  sweeps  across 
the  landscape  and  sends  back  numerous  measurements  with  precise  x,  y,  and  z 
coordinates,  each  having  an  associated  RGB  color  and  intensity  value.  The 
sophisticated  design  of  the  scanner  enables  point  clouds  to  be  captured  that 

1  Named  “Dirt  Site”  by  Army  test  bed  sponsor. 
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correspond  to  true  point  positions  where  the  laser  pulse  hits  the  object.  The  point 
cloud  represents  the  shape  and  position  of  the  objects  scanned  relative  to  the 
position  of  the  scanner  (Leica  Geosystems  HDS,  Inc.  2004).  See  Table  1  for  a 
complete  summary  of  the  scanner  specifications. 


Table  1 

Specifications  for  the  Leica  HDS3000  3D  Laser  Scanner 

Field  of  View 

360  deg  H  x  270  deg  V 

Positional  Accuracy 

6  mm  at  50  m 

Wavelength 

523  nm 

Spot  Size 

<  6  mm  at  50  m 

Pulse  Rate 

1,000  points/sec  effective  to  100  m 

Maximum  Sample  Density  (spacing) 

1.2  mm 

The  operating  software  used  in  conjunction  with  the  laser  scanner  during  data 
acquisition  was  Cyclone  Version  5.0.  The  Cyclone  software  assists  the  surveyor 
in  capturing  point  cloud  data,  visually  interpreting  and  processing  these  data, 
then  integrating  the  collected  information  into  useful  geospatial  formats. 

When  feasible,  the  scanner  was  positioned  at  approximately  one-quarter 
length  intervals  along  each  field’s  long  side  in  an  effort  to  provide  uniform  scan 
coverage  of  the  entire  site.  Also,  to  obtain  an  unobstructed  viewing  angle  from  all 
scanner  locations,  the  scanner  was  positioned  on  a  trailer-mounted,  chain-driven 
lift  and  raised  to  a  scan  height  of  7.62  m  above  the  ground.  At  each  initial  fixed 
scanner  position,  a  high-resolution  picture  image  of  the  viewing  area  was  cap¬ 
tured  using  the  scanner’s  built-in  camera.  This  allowed  the  surveyor  to  easily 
view  and  depict  areas  to  be  scanned,  making  Cyra  target  acquisition  and  site/ 
minefield  scans  much  more  efficient. 


Cyra  target  acquisition 

After  a  suitable  image  was  acquired  at  each  scanner  setup,  all  corresponding 
Cyra  targets  within  the  effective  scan  range  (<  100  m)  were  probed  with  the 
scanner  prior  to  acquisition  to  determine  the  approximate  distance  to  the  target. 
High-density  laser  scans  (range  0.2  to  2.0  cm  spacing)  were  used  for  Cyra  target 
acquisition  and  were  commensurate  with  size  and  distance  to  the  target  from  each 
scanner  setup.  The  tie  points  generated  from  each  target  acquisition  (Figure  1) 
established  a  set  of  constraints  that  were  used  to  register,  or  geometrically  align, 
all  of  the  scanner  locations  into  a  single,  fitted  point  cloud  representing  each  site. 
Each  tie  point  was  labeled  with  a  unique  registration  ID. 

To  obtain  a  tie  point  with  minimal  deviation  from  the  center  of  the  target,  the 
target  must  be  scanned  with  a  sufficient  density  of  postings  in  the  object’s  center. 
Once  a  post  point  is  manually  selected  that  is  close  to  the  target’s  center,  the 
scanner  performs  a  coarse  scan  to  locate  the  center  circle  then  proceeds  with  a 
fine  scan  (~1.2  mm  spacing  at  50  m)  using  an  algorithm  to  locate  the  exact 
center. 
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Figure  1.  Tie  point  generated  from  precise  Cyra  target  acquisition.  Posting  color 
represents  multi-hue  intensity  of  each  laser  return.  High  intensity 
appears  blue  and  low  intensity  appears  red 

A  vertex  is  placed  at  the  perceived  center  of  the  Cyra  target,  representing  the  tie 
point.  A  minimum  of  three  Cyra  targets  was  placed  at  each  scanner  setup  to 
produce  sufficient  tie  points  to  be  used  as  constraints  in  the  registration  process. 
Each  Cyra  target  was  strategically  positioned  at  opposite  extents  of  the  minefield 
site  and  served  as  common  targets  to  additional  scanner  setups  in  order  to  suc¬ 
cessfully  reference  each  point  cloud  together. 


Site/minefield  scans 

After  target  acquisition  was  complete  at  each  scanner  setup,  the  sites  were 
scanned.  The  laser  scanner  was  positioned  at  four  locations  at  the  Grass  Site,  two 
on  the  western-most  side  and  two  on  the  eastern-most  side  of  the  minefield. 
Restricted  access  on  the  north  side  of  the  Dirt  Site  prevented  laser  scan  setups  on 
this  side.  As  a  result,  the  scanner  was  positioned  only  on  the  south  side  and 
provided  six  total  setups  on  the  Dirt  Site.  Moderate-density  scans  were  chosen 
for  point  cloud  generation  of  each  site  with  approximately  10  percent  edge 
overlap  at  each  field  extent  to  capture  the  entire  area  of  interest.  Also,  a  10 
percent  scan  overlap  between  site  scans  ensured  sufficient  point  cloud  data 
collection  for  each  site.  The  Grass  Site  and  the  Dirt  Site  were  scanned  at 
approximate  post  spacings  of  5  and  2  cm,  respectively,  in  the  horizontal  and 
vertical  direction.  Figure  2  depicts  a  small  portion  of  the  scanned  Dirt  Site  and 
illustrates  the  point  cloud  representation  of  the  various  mines  and  fiducial 
markers  in  the  minefield. 
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Figure  2.  Graphical  depiction  of  scanned  Dirt  Site  showing  M20  surface  mines  and  fiducial  markers. 
True-color  laser  postings  are  spaced  every  2  cm 


Data  Analysis 

Once  data  acquisition  was  completed  for  both  sites,  processing  of  the  laser 
data  began.  To  generate  a  three-dimensional,  continuous  point  cloud  representa¬ 
tion  of  each  site,  each  scanner  setup  location’s  point  cloud,  or  ScanWorld,  had  to 
be  coregistered  with  one  another  (interim  registration)  and  referenced  to  a  com¬ 
mon  coordinate  system  (georegistration)  for  use  with  other  spatial  datasets  and  to 
perform  additional  analyses.  A  ScanWorld  can  be  defined  as  a  collection  of 
scanned  point  clouds  that  are  derived  from  consecutive  scans  at  the  same  scanner 
location.  The  ScanWorlds  were  aligned  together  to  form  a  referenced  dataset 
representative  of  a  particular  site.  The  fully  georegistered  scan  data  were  then 
processed  to  produce  digital  elevation  and  terrain  models  applicable  to  each  site. 
Further  descriptions  of  these  post-processing  techniques  are  detailed  below. 
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Laser  data  registration 

Registration  is  a  method  that  aligns  many  individual  ScanWorlds  into  a 
single  georeferenced  ScanWorld  to  represent  the  entire  area  of  interest,  in  this 
case  each  minefield  site.  The  registration  process  makes  use  of  various  mathe¬ 
matical  algorithms  that  compute  the  optimal  alignment  transformations  for  each 
ScanWorld  in  the  registration  such  that  the  constrained  objects  or  point  clouds 
are  aligned  as  closely  as  possible  in  the  resulting  ScanWorld  (Leica  Geosystems 
HDS,  Inc.  2004).  Because  sufficient  Cyra  targets  were  acquired  for  all 
ScanWorlds  at  each  site,  target-based  registration  was  used  for  interim  point 
cloud  alignment.  The  Cyra  targets  placed  at  the  extent  of  each  minefield  served 
as  control  or  tie  points  in  the  registration  process.  These  tie  points  that  were 
common  to  adjacent  ScanWorlds,  appearing  to  be  in  the  same  Cartesian  location 
(x,  y,  and  z  position),  were  fitted  together  to  establish  an  accurate  relationship 
between  each  of  the  ScanWorld  point  clouds.  During  this  initial  registration 
process,  the  Cyclone  software  added  tie  points  as  constraint  objects  to  be  paired 
with  corresponding  tie  points  in  other  ScanWorlds.  The  software  performs  a 
constraint-searching  algorithm  that  locates  objects  with  the  same  registration  ID, 
or  tie  points  that  are  geometrically  consistent,  to  find  the  optimal  solution. 

After  the  ScanWorlds  were  registered  together  for  each  site,  the  single  point 
cloud  was  georegistered  to  the  Universal  Transverse  Mercator  (UTM)  projection 
on  the  basis  of  the  North  American  Datum  1983  (NAD83).  Surveyed  coordinates 
of  specific  Cyra  targets  used  in  the  registration  process  were  collected  to  milli¬ 
meter  accuracies  with  a  Real  Time  Kinematic  (RTK)  Global  Positioning  System 
(GPS)  at  each  site  and  were  used  to  register  the  existing  point  cloud  data  to  the 
UTM  projected  coordinate  system.  Each  target’s  Cartesian  position  was  identi¬ 
fied  to  minimum  accuracies  of  10  mm  horizontal  and  15  mm  vertical.  At  least 
four  surveyed  Cyra  targets  with  known  coordinates  were  used  for  each  site,  and 
these  corresponded  to  all  relative  ScanWorlds.  A  text  file  containing  the 
Northing,  Easting,  and  Elevation  measurements  associated  with  each  surveyed 
Cyra  target  was  imported  into  the  Cyclone  software.  The  imported  survey 
coordinates  served  as  a  new  survey  control  ScanWorld  to  which  all  others  ere 
then  georegistered.  The  surveyed  Cyra  targets  were  positioned  to  achieve  a  large 
aspect  ratio  and  thus  provide  an  optimal  geometric  solution  in  the  registration 
process.  After  successful  georegistration  was  completed  for  both  sites,  surface 
analysis  techniques  were  employed  to  further  characterize  the  sites  and  their 
associated  backgrounds. 


Surface/terrain  analysis 

Topographic  derivatives  were  generated  for  each  site  to  effectively  relate  the 
scanned  laser  data  to  terrain  features.  The  abundance  of  data  points  generated 
from  the  laser  scanner  allowed  for  the  production  of  detailed  digital  terrain  and 
elevation  model  representations  of  each  site.  These  terrain  models  not  only 
provide  3D  visualization  of  the  background  phenomenology  but  also  enable 
analysts  to  measure  topographic  variations  within  the  minefield. 
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The  georegistered  data  points  representing  the  area  of  interest  at  each  site 
were  extracted  and  unified  as  a  single  point  cloud  for  processing.  Scan  data 
outside  the  fence  area,  or  beyond  the  extent  of  each  minefield,  were  discarded 
from  the  data  set  prior  to  processing  to  remove  trees  and  other  superfluous 
background  data.  Tall  vegetation  within  the  Dirt  Site  was  manually  extracted  by 
similar  means  for  generation  of  a  bare  earth  model.  Subsequent  generation  of  a 
contour  map  from  the  bare  earth  model  was  produced  for  the  Dirt  Site  as  well.  In 
addition,  a  vegetation  height  model  was  produced  from  the  laser  point  cloud 
representing  the  Grass  Site. 

The  Cyclone  software  was  used  to  select  five  to  nine  individual  laser  data 
points  representative  of  relatively  flat,  bare  ground  from  a  centralized  area  within 
the  Dirt  Site  point  cloud.  Data  for  these  sites  were  used  as  input  for  a  region 
grow,  surfacing  algorithm.  The  surface-smoothing  algorithm  segments  the  point 
cloud  to  form  a  horizontally  expanded,  planar  point  cloud  indicative  of  the  terrain 
geometry.  The  algorithm  operates  based  on  fit  calculation  parameters  that  are 
user-specified  and  continues  until  all  assumed  non-ground  data  points  are  effec¬ 
tively  isolated  from  the  remaining  ground  points.  The  primary  surface  parameters 
involved  in  this  process  include  (a)  region  thickness  threshold,  which  defines  the 
range  of  data  points  to  be  surfaced  as  ground,  (b)  surface  angle  tolerance  to 
account  for  areas  of  high  relief,  and  (c)  gap  distance,  which  defines  the  maximum 
distance  allowed  between  parts  of  the  same  smooth  surface.  The  region  grow 
algorithm  did  not  properly  identify  certain  points  within  the  Dirt  Site.  Therefore, 
these  point  data  were  manually  edited  until  satisfactory  results  were  obtained. 

After  all  assumed  vegetation  was  removed  and  the  ground  surface  points 
were  identified  for  the  Dirt  Site,  the  points  representing  bare  ground  were  used  to 
create  a  Triangulated  Irregular  Network  (TIN)  or  mesh.  By  producing  a  TIN  of 
the  assumed  ground,  a  coherent  modeled  surface  can  be  easily  visualized.  An 
elevation  contour  map  was  subsequently  produced  from  the  TIN  for  the  Dirt  Site. 
Major  contours  were  specified  at  0.5-m  intervals,  and  the  number  of  minor  inter¬ 
vals  per  major  contour  was  set  at  five.  As  a  result,  this  produced  index  contours 
at  half-meter  intervals  and  a  contour  interval  of  10  cm,  effectively  yielding  a 
highly  detailed,  micro  topographical  profile  of  the  Dirt  Site.  A  regularly  spaced 
sample  grid  was  then  generated  from  the  original  TIN  layer  to  provide  a  digital 
terrain  model  of  the  Dirt  Site. 

Due  to  the  dense  vegetation  present  at  the  Grass  Site,  a  vegetation  height 
map  was  produced  to  better  quantify  the  background  component  of  the  site.  The 
laser  data  points  representative  of  the  Grass  Site  were  exported  from  the  Cyclone 
software  as  an  x,  y,  and  z  text  file.  This  text  file  was  imported  into  a  custom, 
proprietary  application  written  specifically  for  this  vegetation  height  extraction 
(personal  communication,  R.  E.  Melton,  Jr.,  Senior  Programmer,  JAY  A 
Corporation  2004).  The  application  was  designed  to  distinguish  and  isolate 
assumed  ground  hits  and  maximum  vegetation  height  points.  The  application 
extracts  laser  data  from  the  lowest  10  percentile  using  each  point’s  elevation  (z) 
value  and  then  averages  those  within  a  one  square  meter  cell  size.  This  is  the 
assumed  ground.  Likewise,  laser  data  from  the  top  10  percentile  were  extracted 
by  z  value,  averaged,  and  then  output  as  a  single  point  representing  the  average 
maximum  vegetation  height  for  that  1-sq-m  cell.  The  output,  x  and  y  values  for 
the  center  of  each  cell  and  an  average  elevation  value,  were  uploaded  into  an 
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ESRI  point  shapefile.  Vegetation  height  was  calculated  by  subtracting  the 
assumed  ground  elevation  points  from  the  assumed  vegetation  elevation  points. 
These  new  elevation  point  values  were  the  representative  vegetation  height  value 
for  each  one-meter  cell  center,  which  were  used  to  generate  a  1-m  grid  of  the 
entire  Grass  Site. 
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3  Results  and  Discussion 


Laser  Data  Registration 

Mean  Absolute  Error  (MAE)  was  used  to  measure  the  accuracy  of  the  point 
cloud  registrations  for  each  site  (Table  2).  Point  cloud  alignment  error  was 
evaluated  for  both  the  interim  registration  and  georegistration  to  the  UTM 
projected  coordinate  system  using  the  RTK-GPS  collected  survey  coordinates. 


Table  2 

MAE  Values  for  Interim  Point  Cloud  Registration  and 

Georegistration  for  Both  Sites 

Minefield  Site 

Registration  Error  (MAE) 

#  ScanWorlds 

Interim 

Geo (UTM) 

Dirt  Site 

0.008  m 

0.014  m 

6 

Grass  Site 

0.014  m 

0.017  m 

4 

MAE  can  be  defined  as  the  weighted  average  of  the  absolute  errors,  with  the 
relative  frequencies  as  the  weight  factors.  Additionally,  the  minimum  MAE  value 
can  be  interpreted  as  the  mean  absolute  deviation  among  data  points.  All  tie  point 
constraints  were  equally  weighted  during  the  registration  process.  One  constraint 
pair,  a  target  tie  point  (DBM  1)  and  its  corresponding  survey  coordinate,  was 
disabled  in  the  final  georegistration  of  the  Dirt  Site  because  of  the  inordinate 
error  compared  to  all  other  constraints  (Figure  3).  See  Appendix  A  for  complete 
registration  diagnostics  for  both  minefield  sites. 

The  Positional  error  of  the  Leica  scan  data  alone  is  generally  around  6  mm 
(Leica  Geosystems  HDS,  Inc.  2004).  However,  the  overlap  measurements  are 
often  imprecise,  especially  for  scans  of  complex  geometry  such  as  grass  and 
other  vegetation.  Various  reasons  exist  that  could  explain  this  intricacy.  Over¬ 
lapping  laser  points  are  often  not  always  from  the  same  surface,  or  blade  of  grass. 
Multiple  viewpoints  from  different,  surrounding  ScanWorlds  generate  a  some¬ 
what  convoluted  perspective  of  the  same  spatial  area,  particularly  for  grass  or 
vegetation  scans,  therefore  increasing  the  point-to-point  deviation  or  error.  Other 
reasons  may  exist  including  environmental  factors  such  as  wind  or  sun  angle, 
which  may  cause  vegetated  surfaces  to  be  more  spatially  variable  between 
consecutive  scans,  contributing  to  a  higher  MAE.  Moderate  winds  during  a 
portion  of  the  Grass  Site  data  capture  contributed  to  lift  platform  and  scanner 
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Figure  3.  Illustration  of  Cyclone  georegistration  process  for  Dirt  Site  depicting  synchronized  windows  of 
corresponding  ScanWorlds  with  disabled  tie  point  constraint  (DBM1) 


sway,  making  precise  tie  point  acquisition  much  more  difficult.  Also,  radiant  heat 
energy  from  the  intense  mid-day  sun  during  the  same  scanning  operation  pro¬ 
duced  an  evident  “wavy”  pattern  in  a  portion  of  the  scan  data.  It  is  speculated, 
therefore,  that  these  were  the  primary  known  causes  of  the  higher  error  for  the 
Grass  Site  (0.017  m)  when  compared  to  the  Dirt  Site  (0.014  m). 


Surface/Terrain  Analysis 

Terrestrial  3D  laser  characterization  efforts  of  the  variable-surface  sites 
yielded  an  elevation  contour  map  and  vegetation  height  model  of  the  Dirt  Site 
and  Grass  Site,  respectively.  Figure  4  illustrates  the  resulting  ground  surface 
elevation  contour  map,  with  contour  intervals  of  10  cm. 
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Figure  4.  Gound  elevation  contour  map  (10-cm  Interval)  detailing  the  micro  topography  representative  of 
the  Dirt  Site 

The  abundance  of  data  points  generated  from  the  laser  scanner  provided  a 
very  rich  data  set  from  which  to  produce  a  very  detailed  micro  topographical 
contour  map  of  the  Dirt  Site.  By  producing  a  10-cm  contour  interval  map,  the 
vertical  relief  and  landscape  profile  of  the  site  were  accurately  depicted.  The 
contour  interval  measurement  chosen  produced  a  minimum  vertical  distance 
between  adjacent  contour  lines,  allowing  for  precise  modeling  capabilities  of  the 
terrain  surface. 

A  vegetation  height  model,  generated  from  the  closely  spaced  laser  data 
points,  effectively  characterized  the  very  dense  vegetative  component  of  the 
Grass  Site  (Figure  5).  Notice  the  high  tree  crown  tops  represented  in  red  on  the 
extreme  northern  and  southwestern  parts  of  the  Grass  Site.  Inspection  of  the 
resulting  vegetation  surface  revealed  the  variation  in  vegetation  height 
differences  across  the  site. 
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Figure  5.  Vertical  view  of  Grass  Site  depicting  vegetation  height  in  meters.  The  image  portrays  a 
chromatic  sequence  with  lower  vegetation  heights  appearing  blue  and  higher  vegetation 
heights  appearing  red 
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4  Future  Considerations 


A  primary  objective  of  this  research  was  to  characterize  various  backgrounds 
(sand,  rock,  grass,  soil,  etc.)  typically  present  in  a  minefield  and  evaluate  the 
effectiveness  of  utilizing  a  laser  scanning  device  to  accomplish  this  task.  It  may 
be  apparent  with  additional  research  that  a  smaller  area  can  be  scanned  to 
produce  similar  results  and  satisfy  program  objectives  while  minimizing 
registration  error.  If  successful,  this  would  significantly  reduce  the  amount  of 
data  collected  and  save  a  great  amount  of  time.  Registration  error  could  also  be 
reduced  by  restricting  Cyra  target  and  laser  measurement  acquisitions  to  within 
75  m,  well  within  the  stated  effective  range  of  the  scanner  (100  m). 

Due  to  the  oblique,  off-nadir  measurement  angle  of  the  elevated  scanner, 
vertical  measurements  of  scanned  objects  are  not  effectively  obtained.  To  better 
determine  the  terrestrial  scanner’s  ability  to  accurately  extract  vegetation  heights 
in  a  minefield,  a  ground-truth  exercise  should  be  implemented  to  develop  a 
control  measurement  of  vegetation  heights  to  compare  to  the  scanner  data 
estimates.  In  addition,  the  related  estimation  of  assumed  “ground”  at  a  vegetated 
site  using  an  average  of  z  values  in  the  lowest  tenth  percentile  of  points  for  a  unit 
area  is  a  methodology  that  should  be  field  validated. 
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Appendix  A 

Registration  Diagnostics 


status:  valid  Registration 
Mean  Absolute  Error: 

for  Enabled  Constraints  -  0.014  m 
for  Disabled  Constraints  •  0.100 
Date:  2004.09.10  09:00:10 
Database  name  :  Dirt_Site 

Scanworlds 
scanworld  1 
Scanworld  2 
Scanworld  3 
Scanworld  4 
Scanworld  S 
scanworld  6 
BM-xyz-os 


Final  Registration  Di agnos tics_UTMxyzDS . txt 


t5 

Scanworld 

Scanworld 

Type 

On/Off 

DBMl 

Scanworld 

1 

Scanworld 

2 

Coincident 

vertex-vertex 

On 

DBM  2 

Scanworld 

1 

Scanworld 

2 

coincident 

vertex-vertex 

On 

Target2 

Scanworld 

1 

Scanworld 

2 

Coincident 

vertex-vertex 

On 

Target3 

scanworld 

2 

Scanworld 

3 

Coincident 

vertex-vertex 

On 

DBM4 

Scanworld 

2 

Scanworld 

3 

Coincident 

vertex-vertex 

On 

DBM  5 

Scanworld 

2 

Scanworld 

3 

Coincident 

vertex-vertex 

On 

Taraet3 

scanworld 

2 

scanworld 

4 

Coincident 

vertex-vertex 

On 

dbm4 

Scanworld 

2 

Scanworld 

4 

Coincident 

vertex-vertex 

On 

DBMS 

Scanworld 

2 

Scanworld 

4 

Coincident 

vertex-vertex 

On 

Target3 
DBM  4 

scanworld 

3 

Scanworld 

4 

coincident 

vertex-vertex 

On 

Scanworld 

3 

Scanworld 

4 

coincident 

vertex-vertex 

On 

DBMS 

scanworld 

3 

Scanworld 

4 

Coincident 

vertex-vertex 

On 

DBM6a 

scanworld 

4 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

DBM6b 

Scanworld 

4 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

DBM  7 

scanworld 

4 

scanworld 

5 

Coincident 

vertex-vertex 

On 

DBM8 

Scanworld 

4 

scanworld 

5 

Coincident 

vertex-vertex 

On 

DBMSa 

scanworld 

4 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

DBM6b 

Scanworld 

4 

Scanworld 

6 

coincident 

vertex-vertex 

on 

DBM  7 

Scanworld 

4 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

DBM8 

scanworld 

4 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

DBM  7 

Scanworld 

5 

Scanworld 

6 

Coincident 

Vertex-Vertex 

On 

DBM8 

scanworld 

5 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

DBM6a 

scanworld 

5 

scanworld 

6 

coincident 

vertex-vertex 

On 

DBM6b 

Scanworld 

5 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

DBMl 

Scanworld 

1 

BM-xyz-DS 

Coincident 

vertex-vertex 

On 

DBM2 

Scanworld 

1 

BM-xyz-DS 

Coincident 

vertex-vertex 

On 

DBMl 

Scanworld 

2 

BM - xyz  DS 

Coincident 

vertex-vertex 

Off 

DBM2 

scanworld 

2 

BM-xyz-os 

Coincident 

vertex-vertex 

On 

DBM  3  a 

Scanworld 

2 

BM-xyz-DS 

Coincident 

vertex-vertex 

On 

0BM3b 

scanworld 

2 

BM-xyz-DS 

Coincident 

vertex-vertex 

On 

DBM9a 

Scanworld 

2 

BM-xyz-DS 

coincident 

vertex-vertex 

On 

DBM9b  scanworld 

Transformations 

2 

BM-xyz-DS 

Coincident 

vertex-vertex 

On 

Name 

rargetlO 

TargetiD 

TargetiD 

TargetID 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetTD 

TargetiD 

TargetiD 

TargetlO 

Targetio 

TargetiD 

TarqetID 

Targetio 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

TargetiD 

Targetio 

TargetiD 

TargetiD 

Targetio 


Scanworld  1 

translation:  (417S364. 153. 
rotation:  (-0.0043,  0.0047, 

Scanworld  2 

transl at ion :  (4175402.145.  577100.966,  342.303)  m 
rotation:  (-0.0073,  0.0047,  - 1 .0000) : 69 . 085  deg 
Scanworld  3 

translation:  (4175430.434,  577218.713,  340.909)  m 
rotation:  (0.0009,  0.0074,  -1 .0000) : 168 . 915  deg 
scanworld  4 

translation:  (4175431.397,  577219.363,  340.954)  m 
rotation:  (0.0023,  -0.0032,  1 . OOOO) : 77 . 021  deg 
Scanworld  5 

translation:  (4175451. 473 ,  577313.238,  346.698)  m 
rotation:  (-0.0079,  0.0115,  -0.9999) : 92 .921  deg 
Scanworld  6 

translation:  (417S498. 617,  577303.640.  347.103)  m 
rotation:  (-0.0079.  0.0111,  -0.9999) : 122.630  deg 
Scanworld  bm  xyz  os 

translation:  (0.000,  0.000,  0.000)  m 
rotation:  (0.0000.  1.0000.  0 .0000) : 0 . 000  deg 

unused  control  space  objects 
Scanworld  1: 

vertex  :  TargetiD  :  Targetl 
Scanworld  3: 

vertex  :  TargetiD  :  dbm4 
scanworld  4 
vertex 


577001.415,  342.135) 
-1.0000) : 71.050  deg 


weight 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1 . 0000 
1.0000 
1.0000 
1 . 0000 
1.0000 
1.0000 
1.0000 
1 . 0000 
1.0000 
1.0000 
1.0000 
1 . 0000 
1.0000 
1.0000 
1 . 0000 
1.0000 
1.0000 
1.0000 
1 . 0000 
1.0000 
1.0000 
1 . 0000 
1.0000 
1 . 0000 
1 . 0000 
1.0000 


Error  vector 

(0.048.  -0.006,  -0.008)  m 
(-0.008,  -0.009,  0.007)  m 
(-0.015,  0.009,  -0.001)  m 
(0.001,  0.003,  0.000)  m 
(0.001.  -0.005,  0.000)  m 
(-0.002,  0.002,  0.000)  m 
(-0.003,  -0.001,  0.000)  m 
(0.006,  0.000,  0.000)  m 
(-0.003,  0.001,  0.000)  m 
(-0.004,  -0.004,  0.000)  m 
(0.005,  0.005,  0.000)  m 
(-0.001,  -0.001,  0.000)  m 
(-0.002,  0.004.  0.002)  m 
(0.000,  -0.008,  0.006)  m 
(0.002,  0.002,  -0.004)  m 
(0.000,  0.002,  -0.004)  m 
(-0.006,  0.000,  0.001)  m 
(-0.002,  -0.014,  0.005)  m 
(0.007,  0.007,  -0.003)  m 
(0.001,  0.008,  -0.004)  m 
(0.005,  0.005,  0.001)  m 
(0.001,  0.005,  0.001)  m 
(-0.004,  -0.004,  0.000)  m 
(-0.002,  -0.006.  -0.001)  m 
(0.049,  0.011,  6.008)  m 
(0.022,  -0.005,  -0.006)  m 
(0.097,  0.017,  0.016)  m 
(0.031,  0.004,  -0.013)  m 
(-0.049.  -0.014,  0.000)  m 
(-0.049,  -0.012,  0.004)  m 
(-0.002,  0.007,  0.004)  m 
(-0.002,  0.009,  0.004)  m 


TargetiD  :  Target  4 


Figure  A1 .  Final  Registration  Diagnostics_UTMxyzDS.txt 


Appendix  A 


Registration  Diagnostics 


interim  Registration  oiagnostics_DS.txt 


Status:  valid  Registration 
Mean  Absolute  Error: 

for  enabled  Constraints  *  0.008  m 
for  Disabled  Constraints  -  0.000  m 
Date:  2004.08.19  13:50:10 
Database  name  :  Dirt_Site 

Scanworlds 
Scanworld  1 
scanworld  2 
Scanworld  3 
Scanworld  4 
Scanworld  5 
scanworld  6 

Constraints 


DBMl 

scanworld 

Scanworld 

1 

scanworld 

Scanworld 

2 

Type 

Coincident:  vertex-vertex 

On/Off 

On 

weight 

1.0000 

Error 

0.024 

m 

DBM  2 

Scanworld 

1 

Scanworld 

2 

Coincident:  vertex  vertex 

On 

1.0000 

0.017 

m 

Target2 

Scanworld 

1 

Scanworld 

2 

Coincident:  vertex-vertex 

On 

1 . 0000 

0.012 

m 

Target3 

DBM4 

scanworld 

2 

Scanworld 

3 

Coincident:  vertex-vertex 

On 

1.0000 

0.003 

m 

Scanworld 

2 

Scanworld 

3 

Coincident:  vertex-vertex 

On 

1.0000 

0.005 

Nl 

DBMS 

Scanworld 

2 

Scanworld 

3 

Coincident:  vertex-vertex 

On 

1.0000 

0.003 

m 

Target3 

scanworld 

2 

Scanworld 

4 

Coincident:  vertex-vertex 

On 

1.0000 

0.003 

m 

DBM4 

scanworld 

2 

scanworld 

4 

coincident:  vertex-vertex 

On 

1.0000 

0.006 

m 

DBMS 

Scanworld 

2 

Scanworld 

4 

Coincident:  vertex-vertex 

On 

1.0000 

0.003 

m 

Target3 

DBM4 

scanworld 

3 

Scanworld 

4 

Coincident:  vertex-vertex 

On 

1.0000 

0.006 

m 

Scanworld 

3 

Scanworld 

4 

Coincident:  vertex-vertex 

On 

1.0000 

0.007 

m 

DBMS 

scanworld 

3 

Scanworld 

4 

Coincident:  vertex-vertex 

On 

1.0000 

0.001 

m 

DBM6a 

Scanworld 

4 

Scanworld 

5 

Coincident:  vertex-vertex 

On 

1.0000 

0.005 

m 

DBM6b 

Scanworld 

4 

Scanworld 

5 

Coincident:  vertex-vertex 

On 

1.0000 

0.010 

m 

DBM7 

Scanworld 

4 

scanworld 

5 

Coincident:  vertex-vertex 

on 

1.0000 

0.005 

m 

DBM8 

scanworld 

4 

scanworld 

5 

Coincident:  vertex-vertex 

On 

1.0000 

0.005 

m 

DBM6a 

Scanworld 

4 

Scanworld 

6 

Coincident:  vertex -vertex 

On 

1.0000 

0.006 

m 

DBM6b 

scanworld 

4 

Scanworld 

6 

Coincident:  vertex-vertex 

On 

1 . 0000 

0.016 

m 

DBM7 

Scanworld 

4 

Scanworld 

6 

Coincident:  vertex-vertex 

On 

1.0000 

0.010 

m 

DBMS 

scanworld 

4 

Scanworld 

6 

Coincident:  vertex-vertex 

On 

1.0000 

0.009 

m 

DBM  7 

Scanworld 

5 

Scanworld 

6 

Coincident:  vertex-vertex 

On 

1.0000 

0.007 

m 

D8M8 

scanworld 

S 

Scanworld 

6 

Coincident:  vertex-vertex 

On 

1.0000 

0.006 

m 

DBM6a 

Scanworld 

5 

Scanworld 

6 

Coincident:  vertex-vertex 

On 

1.0000 

0.006 

m 

DBM6b 

scanworld 

5 

Scanworld 

6 

Coincident:  vertex-vertex 

On 

1.0000 

0.007 

m 

TargetlD: 

TargetID: 

TargetlD: 

TargetIO: 

TargetlD: 

TargetIO: 

TargetlD: 

TargetlD: 

TargetlD: 

TargetIO: 

TargetlD: 

TargetlD: 

TargetIO: 

TargetIO: 

TargetlD: 

TargetlD: 

TargetlD: 

TargetlD: 

TargetIO: 

TargetlD: 

TargetlD: 

TargetlD: 

TargetlD: 

TargetlD: 

Scanworld  Transformations 
Scanworld  1 

translation:  (0.000,  0.000,  0.000)  m 
rotation:  (0.0000,  1.0000,  0 . 0000) : 0 . 000  deg 
Scanworld  2 

translation:  (-3.629,  106.503,  0.952)  m 
rotation:  (-0.0673,  -0.0666,  0.9955) :1. 980  deg 
Scanworld  3 

translation:  (7.843,  227.063,  0.485)  m 
rotation:  (-0.0007.  0.0117,  -0.9999) : 97. 860  deg 
Scanworld  4 

translation:  (7.142,  227.990,  0.53S)  m 
rotation:  (0.0027,  -0.0049,  1. 0000) : 148 . 077  deg 
Scanworld  5 

translation:  (18.621.  323.252.  7.015)  m 
rotation:  (-0.0297,  6.0179,  -0.9994) : 21. 872  deg 
Scanworld  6 

translation:  (-29.086,  329.476,  7.388)  m 
rotation:  (-0.0173,  0.0111.  -0.9998) : 51. 577  deg 

unused  Control  Space  Objects 
Scanworld  1: 

vertex  :  TargetlD  :  Target 1 

Scanworld  2: 
vertex  : 
vertex  : 
vertex  : 
vertex  : 


TargetlD  :  DBM3a 
TargetlD  :  D8M3b 
TargetIO  :  DBM9a 
TargetlD  :  DBM9b 


Error  vector 

(0.024.  -0.004,  -0.002)  m 
(-0.016,  -0.005,  0.001)  m 
(-0.007,  0.009,  0.001)  m 
(0.000.  0.003,  0.000)  m 
(-0.003,  -0.004,  0.000)  m 
(0.003,  0.001,  0.000)  m 
(0.002.  -0.002,  0.000)  m 
(-0.006,  0.002,  0.000)  m 
(0.003,  0.000,  0.000)  m 
(0.002.  -0.005,  0.000)  m 
(-0.003,  0.006.  0.000)  m 
(0.001.  -0.001,  0.000)  m 
(0.003,  0.003,  0.002)  m 
(-0.003,  -0.008,  0.006)  m 
(-0.001,  0.002,  -0.004)  m 
(0.001.  0.002,  -0.004)  m 
(0.006,  -0.002,  0.001)  m 
(-0.003,  -0.014,  0.005)  m 
(-0.004,  0.009,  -0.003)  m 
(0.001,  0.008,  -0.004)  m 
(-0.003,  0.006,  0.001)  m 
(0.001,  0.006,  0.001)  m 
(0.002,  -0.005,  0.000)  m 
(0.000,  -0.007,  -0.001)  m 


scanworld  3: 
Vertex  : 

Scanworld  4: 
vertex  : 


TargetIO  :  dbm4 
TargetlD  :  Target  4 


Figure  A2.  Interim  Registration  Diagnostics_DS.txt 
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Appendix  A  Registration  Diagnostics 


Final  Registration  oiagnostics_UTMxyzGS.txt 


status:  valid  Registration 
Mean  Absolute  Error: 

for  Enabled  constraints  «  0.017  m 
for  oisabled  constraints  «  0.000 
oate:  2004.09.10  09:11:26 
Database  name  :  Grassy_Sl te_Tal 1 veg 

scanworlds 
scanworld  2 
Scanworld  3 
Scanworld  5 
Scanworld  6 
BM-xyz-GS 

constraints 
Name 

Targetio: 

TargetlD: 

TargetID: 
largetio: 

TargetlD : 

TargetID: 

TargetlD: 

TargetID: 

TargetID: 

TargetID: 

TargetlD: 

Targetio: 

TargetlD: 

Targetio: 

TargetID: 

Targetio: 

Targetio: 

Targetio: 

Targetio: 

TargetID: 

TargetID: 

TargetID: 

Targetio: 

TargetlD: 

Targetio: 

TargetlD: 

TargetID: 

Targetio: 

TargetID: 

Targetio: 

Targetio: 

Targetio: 

Targetio: 

Targetio: 

TargetID: 

Targetio : 

TargetID: 

TargetID: 

TargetlD: 

TargetID: 

TargetID: 

TargetID: 

scanworld  Transformations 
Scanworld  2 

translation:  (4176111.191, 
rotation:  (-0.0025,  0.0015, 

Scanworld  3 

translation:  (4176069.258, 
rotation:  (-0.0009,  0.0019, 
scanworld  5 

translation:  (4176128.055,  576799.248,  343.859)  m 
rotation:  (-0.0089,  -0.0070,  0 . 9999) : 64 . 915  deg 
Scanworld  6 

translation:  (4176170.164,  576836.374,  34S.654)  m 
rotation:  (-0.0114,  -0.0030,  0.99991:54.224  deg 
Scanworld  BM-xyz-GS 
translation:  (0.000,  0.000,  0.000)  m 
rotation:  (0.0000,  1.0000,  0 .0000) : 0 . 000  deg 
Unused  ControlSpace  Objects  :  none 


BMl 

Scanworld 

Scanworld 

2 

Scanwor 1 d 
scanworld 

3 

Type 

Coincident 

vertex-vertex 

on/off 

On 

weight 

1.0000 

Error 

0.022 

Error  vector 
m  (-0.016,  0.010, 

BM2 

scanworld 

2 

scanworld 

3 

Coincident 

vertex-vertex 

on 

1.0000 

0.025 

m  (0.011,  -0.013, 
m  (-0.001,  0.011, 

BM3 

Scanworld 

2 

Scanworld 

3 

Coincident 

vertex-vertex 

On 

1.0000 

0.012 

TargS 

scanworld 

2 

scanworld 

3 

coincident 

vertex-vertex 

On 

1 . oooo 

0.008 

m  (0.006,  -0.003, 

BMl 

Scanworld 

2 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.008 

m  (-0.006,  0.000. 

BM2 

Scanworld 

2 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.005 

m  (-0.003,  -0.001 

BM3 

Scanworld 

2 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.015 

m  (0.005,  0.014,  ( 

Targ5 

Scanworld 

2 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.011 

m  (0.003.  -0.007, 

BMl 

Scanworld 

2 

scanworld 

6 

coincident 

vertex-vertex 

on 

1.0000 

0.007 

m  (0.003,  0.002, 

BM2 

Scanworld 

2 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.009 

m  (-0.006,  -0.002 

BM3 

Scanworld 

2 

Scanworld 

$ 

Coincident 

vertex-vertex 

On 

1.0000 

0.007 

m  (-0.007.  -0.004 

TargS 

Scanworld 

2 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.014 

m  (0.013.  0.004,  ( 

BM3 

Scanworld 

3 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.010 

m  (0.007.  0.002,  ( 

BM2 

Scanworld 

3 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.023 

m  (-0.014.  0.012, 

BMl 

Scanworld 

3 

scanworld 

5 

Coincident 

vertex-vertex 

On 

1 . 0000 

0.016 

m  (0.010,  -0.010, 

BM4 

Scanworld 

3 

scanworld 

5 

coincident 

vertex-vertex 

On 

1 . 0000 

0.005 

m  (0.001,  -0.001. 
m  (-0.003.  -0.004 

TargS 

scanworld 

3 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.013 

BM3 

scanworld 

3 

Scanworld 

6 

Coincident 

Vertex-vertex 

On 

1.0000 

0.016 

m  (-0.005,  -0.015 

BMc 

Scanworld 

3 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.032 

m  (-0.017.  0.011, 

BMl 

Scanworld 

3 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.027 

m  (0.018,  -0.009, 

BM4 

Scanworld 

3 

scanworld 

6 

coincident 

vertex-vertex 

On 

1.0000 

0.015 

m  (-0.002.  0.008, 

TargS 

Scanworld 

3 

Scanworld 

6 

Coincident 

vertex-vertex 

on 

1.0000 

0.010 

m  (0.007,  0.007, 

BM3 

scanworld 

5 

Scanworld 

6 

coincident 

vertex-vertex 

On 

1.0000 

0.021 

m  (-0.012,  -0.017 

BM2 

Scanworld 

5 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.011 

m  (-0.003,  -0.001 

BMl 

scanworld 

5 

scanworld 

6 

coincident 

vertex-vertex 

On 

1.0000 

0.014 

m  (0.009,  0.002, 
m  (-0.003,  0.009. 

BM4 

Scanworld 

5 

Scanworld 

0 

Coincident 

vertex-vertex 

On 

1.0000 

0.013 

TargS 

Scanworld 

S 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.019 

m  (0.010,  0.012,  < 

BMl 

Scanworld 

2 

BM-xyz-GS 

coincident 

vertex-vertex 

On 

1.0000 

0.024 

m  (-0.017.  -0.013 

BM2 

Scanworld 

2 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.021 

m  (0.012,  -0.010, 

BM3 

scanworld 

2 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.011 

m  (0.002,  0.011,  ( 

BM3 

Scanworld 

3 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.009 

m  (0.004,  -0.001, 

8M2 

scanworld 

3 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.032 

m  (0.000,  0.003,  < 
n  (-0.001,  -0.023 

BMl 

Scanworld 

3 

BM-xyz-GS 

coincident 

vertex-vertex 

On 

1.0000 

0.033 

BM4 

scanworld 

3 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.030 

n  (-0.003,  0.025, 

BM3 

scanworld 

5 

BM-xyz-GS 

coincident 

vertex-vertex 

On 

1.0000 

0.004 

m  (-0.003,  -0.003 

BM2 

Scanworld 

5 

BM-xyz-GS 

Coincident 

vertex  vertex 

on 

1.0000 

0.025 

m  (0.015,  -0.009, 
m  (-0.011.  -0.013 

BMl 

scanworld 

5 

BM-xyz-GS 

coincident 

vertex-vertex 

On 

1.0000 

0.023 

BM4 

Scanworld 

5 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.028 

m  (-0.004,  0.026, 

BMl 

Scanworld 

6 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.017 

m  (0.009,  0.014,  1 

BM2 

Scanworld 

6 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.021 

m  (0.018.  -0.009. 

BMl 

Scanworld 

6 

BM-xyz-GS 

Coincident 

vertex-vertex 

On 

1.0000 

0.025 

in  (-0.020,  -0.015 

BM4 

Scanworld 

6 

BM-xyz-GS 

coincident 

vertex-vertex 

On 

1.0000 

0.017 

m  (-0.001.  0.017. 

0.012)  m 
-0.018)  ir 

-0.005)  m 
0.00S)  m 
0.005)  m 
,  -0.004) 
0.002)  m 
-0.007)  nr 
-0.006)  m 
.  0.007)  rr 


0.014)  m 
-0.008)  m 
-0.004)  m 


-0.018)  m 
-0.013)  m 


0.032)  m 
,  -0.024) 
-0.015)  n 
,  0.001)  m 
0.018)  m 
,  -0.016) 
-0.011)  m 
0.003)  m 


-0.002)  n 


576907.227,  351.818)  m 
-1.0000):135.445  deg 

576873.798,  349.791)  m 
-1.0000:130.307  deg 


Figure  A3.  Final  Registration  Diagnostics_UTMxyzGS.txt 


Appendix  A  Registration  Diagnostics 


interim  Registration  Diagnostics_GS.txt 


status:  valid  Registration 

Mean  Absolute  Error: 

for  Enabled  Constraints  -  0.014  m 
for  Disabled  Constraints  =  0.000 
Date:  2004.08.19  12:04:00 
Database  name  :  Grassy_Si te_Tal lveg 

Scanworlds 
scanworld  2 
Scanworld  3 
Scanworld  5 
Scanworld  6 


15 

Scanworld 

Scanworld 

Type 

On/Off 

weight 

1.0000 

Error 

BMl 

Scanworld 

2 

Scanworld 

3 

coincident 

vertex-vertex 

On 

0.020 

m 

BM2 

scanworld 

2 

Scanworld 

3 

Coincident 

vertex-vertex 

on 

1 . OOOO 

0.024 

m 

bm3 

Scanworld 

2 

Scanworld 

3 

Coincident 

vertex-vertex 

On 

1.0000 

0.012 

m 

Targ5 

Scanworld 

2 

Scanworld 

3 

Coincident 

vertex-vertex 

On 

1.0000 

0.010 

m 

BMl 

Scanworld 

2 

scanworld 

5 

coincident 

vertex-vertex 

On 

1.0000 

0.006 

m 

BM2 

Scanworld 

2 

Scanworld 

5 

Coincident 

vertex -vertex 

On 

1.0000 

0.004 

m 

BM3 

scanworld 

2 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1 . 0000 

0.014 

m 

Targ5 

Scanworld 

2 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1 . OOOO 

0.010 

m 

BMl 

Scanworld 

2 

scanworld 

6 

coincident 

vertex-vertex 

On 

1 . 0000 

0.009 

m 

BM2 

scanworld 

2 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1 . 0000 

0.012 

m 

bm3 

Scanworld 

2 

scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.009 

m 

TargS 

scanworld 

2 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.014 

m 

BM3 

scanworld 

3 

scanworld 

5 

coincident 

vertex-vertex 

On 

1 . OOOO 

0.011 

m 

BM2 

Scanworld 

3 

Scanworld 

5 

Coinci dent 

vertex-vertex 

On 

1.0000 

0.023 

m 

BMl 

scanworld 

3 

scanworld 

5 

Coinci dent 

vertex-vertex 

On 

1 . OOOO 

0.016 

m 

BM4 

scanworld 

3 

Scanworld 

5 

Coincident 

vertex-vertex 

On 

1 . 0000 

0.004 

m 

Targ5 

scanworld 

3 

scanworld 

5 

Coincident 

vertex-vertex 

On 

1.0000 

0.013 

m 

BM3 

Scanworld 

3 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.017 

m 

bm2 

scanworld 

3 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1 . OOOO 

0.032 

m 

BMl 

Scanworld 

3 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.027 

m 

BM4 

scanworld 

3 

scanworld 

6 

coincident 

vertex-vertex 

On 

1 . 0000 

0.016 

m 

TargS 

Scanworld 

3 

Scanworld 

6 

Coincident 

Coincident 

vertex-vertex 

On 

1 . 0000 

0.010 

m 

bm3 

scanworld 

5 

Scanworld 

6 

vertex-vertex 

On 

1 . OOOO 

0.022 

m 

bm2 

Scanworld 

5 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.011 

m 

BMl 

scanworld 

5 

scanworld 

6 

Coincident 

vertex-vertex 

On 

1 . 0000 

0.013 

m 

BM4 

Scanworld 

5 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1 . 0000 

0.013 

m 

Targ5  Scanworld 

Transformations 

S 

Scanworld 

6 

Coincident 

vertex-vertex 

On 

1.0000 

0.018 

m 

Name 

TargetlD 

TargetID 

TargetlD 

TargetID 

TargetlD 

TargetID 

TargetID 

TargetID 

TargetID 

TargetID 

TargetID 

TargetlD 

TargetID 

TargetID 

TargetID 

TargetID 

TargetID 

TargetlD 

TargetlD 

TargetID 

TargetlD 

TargetID 

TargetID 

TargetID 

TargetID 

TargetID 

TargetlD 


Scanworld  2 

translation:  (0.000,  0.000, 
rotation:  (0.0000,  1.0000, 

Scanworld  3 
translation:  (53.235,  6.436,  -2.173)  m 
rotation:  (0.0052,  0.0313,  0.9995) :5. 141  deg 
Scanworld  5 

translation:  (65.086.  -87.740,  -8.533)  m 
rotation:  (-0.0030,  6.0083,  -1.0000) : 159. 644  deg 
Scanworld  6 

translation:  (9.095,  -91.707,  -6.572)  m 
rotation:  (-0.0008,  0.0076.  -1. 0000) : 170 . 336  deg 

unused  ControlSpace  Objects  :  none 


0.000)  m 

0.0000) : 0.000  deg 


Error  vector 
(0.004,  0.017,  0.009)  m 
(0.002,  -0.018,  -0.015)  m 
(-0.006,  0.009,  -0.006)  m 
(-0.001,  -0.007,  0.007)  m 
(0.004,  0.003.  0.002)  m 
(0.004,  0.001,  -0.001)  m 
(-0.013,  0.006,  0.003)  m 
(0.004,  -0.007.  -0.004)  m 
(-0.002,  -0.002,  -0.009)  m 
(0.007,  0.002,  0.010)  m 
(0.008,  0.002,  -0.001)  m 
(-0.011,  -0.006,  0.006)  m 
(-0.007,  -0.003,  0.008)  m 
(0.001,  0.019,  0.014)  m 
(0.000,  -0.014,  -0.008)  m 
(0.000,  -0.001,  -0.004)  m 
(0.005.  -0.001,  -0.012)  m 
(0.015,  -0.007,  0.005)  m 
(0.005,  0.020,  0.025)  m 
(-0.006,  -0.019,  -0.018)  m 
(-0.004.  0.007,  -0.014)  m 
(-0.010,  0.000,  -0.002)  m 
(0.021,  -0.004,  -0.003)  m 
(0.003,  0.001.  0.011)  m 
(-0.007,  -0.005.  -0.010)  m 
(-0.004,  0.008,  -0.010)  m 
(-0.015,  0.001,  0.010)  m 


Figure  A4.  Interim  Registration  Diagnostics_GS.txt 
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Appendix  C 

Scanner  Setup  Pictures 


Figure  Cl .  Trailer-mounted  setup  of  Leica  HDS3000  3D  Laser  Scanner  elevated  7.62  m  above  Grass  Site 
(viewed  from  southwest  looking  northeast,  taken  26  July  2004) 
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